Vesicles in a Poiseuille flow 
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Vesicle dynamics in unbounded Poiseuille flow is analyzed using a small-deformation theory. Our 
analytical results quantitatively describe vesicle migration and provide new physical insights. At 
low ratio between the inner and outer viscosity A (i.e. in the tank-treading regime), the vesicle 
always migrates towards the flow centerline, unlike other soft particles such as drops. Above a 
critical A, vesicle tumbles and cross-stream migration vanishes. A novel feature is predicted, namely 
the coexistence of two types of nonequilibrium configurations at the centreline, a bullet-like and a 
parachute-like shapes. 
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The motions of vesicles made of closed bilayer mem- 
branes in viscous flows have received increasing atten- 
tion in recent years because of their relevance to un- 
derstanding cell dynamics in the microcirculation. The 
unusual mechanical properties of the lipid bilayer mem- 
brane, such as fluidity, incompressibihty and resistance to 
bending, give rise to a number of fascinating nonequilib- 
rium features of vesicle micro-hydrodynamics For exam- 
ple, in linear shear flow vesicles are predicted to exhibit 
different types of motion: (i) tank-treading (TT) (the 
fluid membrane rotates as a tank-tread, while the orien- 
tation angle of the vesicle remains fixed in time) P, Q , 
(ii) tumbling (TB) [i,|3], (iii) vaciUating-brcathing (VB) 
(the long axis undergoes oscillation about the flow, while 
the shape shows breathing) Q. Some of these behav- 
iors have been confirmed experimentally d, 0, 111 ■ Vesicle 
dynamics in linear flows is now fairly well understood 

UMMMMM- 

Quadratic flows such as the Poiseuille flow are quite 
common, especially in the blood circulatory system, yet 
the impact of flow curvature on vesicle dynamics is not 
well understood. Number of studies have focused on cap- 
illary flows, where the vesicle diameter is comparable to 
the channel size and lubrication effects control vesicle mo- 
tion [3, [31 1 or vesicles near walls, where hydrodynamic 
interactions with the wall give rise to cross-stream migra- 
tion (in a direction perpendicular to the flow) [l^, 13 ■ 
Blood vessels such as arterioles, however, can be 10-50 
times larger than the cell diameter. In this case, effects 
of flow curvature may become dominant. Vesicle dynam- 
ics in unbounded Poiseuille fiow has been considered only 
to a limited extent, using two-dimensional numerical sim- 
ulations Cross-stream migration has been observed, 
even in the absence of wall [TtI . The physical mechanism 
behind this migration is related to the spatial variations 
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of shear rate that exists in quadratic flows: deformable 
particles tend to move towards regions with lower shear 
in order to minimize shape distortion p^ . Drops, for 
example, migrate towards or away the flow centerline de- 
pending on the viscosity ratio In the case of vesicles, 
however, details of the migration mechanism remain elu- 
sive and there are number of open questions: What is 
the role of flow curvature in cross-stream vesicle migra- 
tion? Is there relation between the migration direction 
and the basic modes of vesicle dynamics (i.e. TT, VB 
and TB)? What physical parameters control vesicle mi- 
gration? The development of a theory that can answer 
these questions represents a challenging task because not 
only the shape, but also the location of the vesicle is not 
known a priori, and it must be solved for in a consistent 
manner. 

Problem formulation and solution: In a reference frame 
centered in the vesicle, a channel fiow is written as 

vo = (-I's - 7y - ay^)e3;, (1) 

where a is a measure of the curvature of the flow profile, 
7 is the local shear rate, which depends to the distance 
between the vesicle center and the flow axis, 7 = 2yoa, 
and Vs is the slip velocity (to be determined), which is 
the difference between the actual velocity of the vesicle 
center in the fiow direction and the unperturbed fiow. 

At the length scales of the vesicle, water is effectively 
very viscous. Hence, the total fiow field v = vq + u (im- 
posed fiow Vo plus perturbation u due to the presence of 
the vesicle) obeys the Stokes equations inside and outside 
the vesicle 

Vp = ryiV^v, V • V = 0, (2) 

where rji [i = 1,2) is the fiuid viscosity; 771,772 are the 
internal and external viscosities, and A = 771/772 is a mea- 
sure of the viscosity contrast. 

Assuming a nearly spherical shape, the solution to 
the creeping-fiow equations is obtained as a regular per- 
turbation expansion in the excess area, which is the 
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difference in the areas of the vesicle and an equiva- 
lent sphere (with radius vq) with the same volume, i.e., 
A = A/vq — An. From the definition of the area it is evi- 
dent that A is a quadratic function of the shape deviation 
from sphere. Thus it is convenient to set e = A^/^ as a 
formal expansion parameter. A time-dependent vesicle 
configuration is described by 



R = Ro + ro[l + e/(t, 



(3) 



where Rg — (xq, yoi -^^o) is the time-dependent position of 
the centroid of the vesicle and / describes shape deviation 
from the spherical shape. 6 and <j) arc spherical coordi- 
nates. The vesicle shape 0) can be decomposed in an 
infinite series of spherical harmonics YimiO,(l)) excluding 
the / = 1 mode (which describes the translation Rq): 
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The amplitudes Fim {t) describe the shape evolution. For 
the present study (leading order analysis) it is sufficient 
to include the I ~ mode, which ensures volume conser- 
vation, and the / = 2, 3 modes, which are excited by the 
imposed Poiseuille flow. 

The velocity field is given by the classical Lamb solu- 
tion [l9l | . This solution contains integration factors which 
are determined from the boundary conditions. These are: 
(i) Continuity of the fiuid velocity across the membrane 
(valid if impermeability is assumed), (ii) Continuity of 
the stress across the membrane: the jump of fluid stress 
across the membrane is balanced by the membrane force. 
The latter consists of a normal force due to resistance to 
bending, and a tangential (fictitious) force, which origi- 
nates from a Lagrange multiplier enforcing local mem- 
brane incompressibility. (iii) Local membrane incom- 
pressibility, which restricts the surface velocity field to be 
solenoidal. We have solved the full hydrodynamic prob- 
lem in the same spirit as in the calculation for the linear 
shear flow. Details can be found in [l^]; here we focus 
solely on the results and discuss the physical implications. 

It should be noted that in certain cases, e.g., drops, 
the migration velocity can be obtained using a simplified 
approach based on the Lorentz reciprocal theorem(20|. 
However, this approach is inapplicable to vesicles because 
of the tangential membrane force. 

Results and Discussion: The time evolution of the 
shape function (dtf) together with the motion of the 
vesicle center with respect to the flow (Ro) are com- 
puted from the kinematic condition, which requires that 
the membrane moves with the fluid velocity. We find for 
the evolution of the 1 — 2,3 (which are the only excited 
modes at leading order) the following equations: 
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with Vt = dt + imayo, and 7V20 = iV2,±i = 0, N2^±2 = 
±8x/30^, 7V30 = iV3,±2 = 0, A^3,±i = T5V2l^/3, 
^^,,±3 ~ T5\/357r. K is the membrane bending rigidity 
rescaled by r-g/??. The evolution equations ^ [H]) con- 
tain the tension-like quantity ctoj which is the homoge- 
neous part of the Lagrange multiplier 0, [l^l . Its value 
is computed from the condition that the shape deforma- 
tion complies with the available excess area, ctq is easily 
determined from the above equations, and the expres- 
sion relating excess area and the shape modes Fim, as 
explained in [l3 |. 

Knowledge of the shape evolution allows us to compute 
the lateral (cross-stream) migration along the ?/-direction 
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This equation is coupled back to the shape, leading to 
some interesting dynamics. First, in contrast to drops, a 
vesicle always migrates towards the centerlinc. Second, 
the migration velocity depends in a non-trivial way on 
the viscosity ratio and excess area. 

For A < Ac, a numerical solution of the coupled equa- 
tions for the shape ^ [S]) and the migration ([7]) yields 
the usual tank-treading ellipsoid as long as the vesicle 
is far from the centerline, because the linear shear dom- 
inates over the quadratic component of the flow. The 
latter, however, induces migration towards the center of 
the flow. As the vesicle moves its shape deforms more and 
more until it becomes a parachute (or assumes another 
shape, namely bullet-like, as discussed below) when the 
centerline is reached (Fig.[T]). The migration velocity far 




FIG. 1: Evolution of the vertical position of the vesicle as a 
function of time. Some snapshots of the vesicle shape along 
the trajectory are shown as insets. Parameters: a = 6, k = 
1, A = 0.25. 

from the centerline is well approximated by plugging the 
solution for the tank-treading ellipsoidal shape, R — e/2, 
cos 2?/' = e/2/i with R and -0 defined by 6^22 = Re'^^'f') 
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into Eq. (0) 1]. 
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For A ~ 1 the migration velocity is approximately Vy = 
— 0.16 aA^/^rg. In a typical experiment in a rectangu- 
lar microcliannel we have a ^ 0.3(//m-s)"\ A = 0.25 
tq = 20 /zm, and a migration velocity of 0.75/xm/s[21 



Using Eq. we find an expected migration velocity of 
about 1 /im/s. Note also that ^ is independent of the 
membrane bending rigidity. This is understood as fol- 
lows. The local shear rate is proportional to the distance 
from the centerline: 7 ~ yo- The elongational component 
of the flow (which increases like the shear rate) must be 
compensated by the membrane tension in order to ful- 
fil local membrane incompressibility. It follows that the 
tension ctq must also scale like yo- Then, for large yo, the 
terms 6k and 12k in Eqs. ([5l [6]) become negligible and 
the membrane dynamics becomes independent of k. 

The direction of the lateral migration of the vesicle in 
Poiseuille flow can intuitively be understood from the fol- 
lowing argument. Assume that we have solved the prob- 
lem for the pure shear flow. 7(2/0)2/- The solution is the 
tank-treading vesicle shown in Fig. Next, we add the 




FIG. 2: (Color online) A vesicle in Poiseuille flow far from the 
centerline. The local velocity field can be decomposed into 
a local shear j{yo) (responsible for the tank-treading vesicle 
at the shown orientation) and a (small) quadratic correction 
(blue arrows). Due to the inclination of the vesicle long axis 
with respect to the quadratic flow field, there is a net force 
acting on the vesicle membrane in the negative y direction. 

quadratic correction —ay^ as a perturbation. The addi- 
tional flow, acting on the membrane of the vesicle, can 
be decomposed into two components: one tangential to 
the membrane, which locally modifies the tank-treading 
velocity, and one normal to the membrane, which mod- 
ifies the vesicle shape and possibly cause migration. In 
the segments AB and CD the perpendicular component 
tends to push the vesicle into the negative y direction, 
whereas in the smaller segments DA and BC, the per- 
pendicular component pushes the vesicle in the positive y 
direction. As the segments AB and CD account for most 
of the vesicle's surface, the net force has a component in 
the negative y direction (towards the flow centerline). 

Let us turn now to the problem of migration in the 
tumbling regime. From the analytical results discussed 
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FIG. 3: Migration velocity as a function of the viscosity con- 
trast A for the various regimes: tank-treading (TT), tumbling 

(TB), anrl vnfil1n + iTnr_V>roafViino- /"VR^ 




FIG. 4: Migration velocity as a function of the excess area 
far from the centerline. Parameters are chosen such that the 
tank-treading regime stops at A = 0.9. 



above for the tank-treading regime, one finds that the 
migration velocity close to the tumbling regime (which 
occurs at A = Ac Q) behaves as — Ac, and thus it 
vanishes exactly at the threshold. This analysis is valid 
on the tank-treading side. For A > Ac, the vesicle tum- 
bles, and in this case we have integrated numerically the 
set of equations (|5l6l7p . The migration (averaged over 
time) is again always towards the flow centerline but the 
magnitude of the migration velocity is very close to zero 
(of the order of 10"^ rg per tumbling cycle). The results 
are shown on FigO For large viscosity ratios, A 1, 
the migration velocity must go to zero, since the vesicle 
behaves as a rigid sphere. The kinematic reversibility of 
the Stokes equations precludes a lateral migration in this 
case. This feature is also evident from the expression for 
the migration velocity ([7]). Indeed, the F/^ s are finite due 
to the constraint of fixed excess area, but the prefactors 
0{X~^). Thus, when A 00, 0. 

An interesting feature is the non-monotonic depen- 
dence of the migration velocity on the excess area A, 
see Fig. 2] For a spherical object (i.e., A = 0) one 
has no migration, = 0, due to the up-down symme- 
try. If viscosity is large enough so that tumbling be- 
comes possible for a sufficiently deflated vesicle, then we 
similarly expect 0, for high enough A. Hence, 
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FIG. 5: Critical angle ^/)c separating the basins of attraction 
for two coexisting solutions at the centerline. 

it follows that the absolute value of the migration ve- 
locity attains a maximum for a certain value of A as 
illustrated in Fig. U) This results also agrees with the 
argument given on the lift force under a linear shear 
fiow[l5|. Vesicle approach to the flow centerline also 
can proceed in a non-monotonous manner, the ratio 
rjoRQa/n (this is a measure of the relative strength of 
hydrodynamic and bending stresses) is sufficiently large 
(typically above 10), the vesicle approaches the centerline 
monotonously. However, for smaller curvature of the flow 
field (or higher membrane rigidity) the vesicle trajectory 
exhibits damped oscillations about the centerline. For 
typical values a ^ 0.1(/ims)~^, Rq ~ lO/im, and by us- 
ing the viscosity of water and k ~ AOksT, one finds that 
rjoRQa/K ^ 0.1. This means these damped oscillations 
should be easily observed experimentally. 

A new surprising feature we have discovered is the co- 
existence of shape solutions at the centerline for small 



enough curvature of the flow field (see Fig. [5]). More 
precisely, this occurs if jyo-Rga/K is lower than about 
4 — 5. For one solution the longest axis of the vesicle 
is oriented in the flow direction (bullet-like shape), and 
for the other one it is perpendicular to it (parachute- 
like shape). Each of these solutions has its own basin 
of attraction, i.e., if the initial angle at the centerline 
is smaller than ipc, the bullet-like shape is attained, and 
otherwise the parachute-like shape. If, however, the cur- 
vature of the flow is strong enough, the vesicle will as- 
sume a pronounced parachute-like shape, and the two 
solutions will be indistinguishable. We hope to report 
further on this matter in a future work. The experimen- 
tal range that is relatively easily accessible is about [lH 
10~^ < TjoRf^a/K < 10. This means that in principle 
this prediction is not devoid of experimental testability. 
The challenge is to orient the vesicle in the appropriate 
basin of attraction. Optical tweezers constitute a possible 
tool for this task. The slip velocity for the parachute-like 
shape is higher, while for a bullet-like vesicle is lower 
than that of a rigid sphere, Vg = — l/3arg (as deduced 
from Faxen's law [l^l)- Since the shp velocity for the 
bullet-shape is smaller than that of the parachute the 
dissipation implied by this morphology should be lower. 

In summary, the dynamical behavior of a vesicle in un- 
bounded quadratic flow has revealed number of novel fea- 
tures that stem from the subtle interplay between mem- 
brane mechanics and shear gradients due to flow cur- 
vature. The analysis has been also performed for an ax- 
isymmetric Poiseuille flow and yields the same qualitative 
results. 
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